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X .  Introduction 

The  problem  we  consider  in  this  paper  is  depicted  in  Figure  1.1.  He 
have  an  abject  that  traverses  a  specified  trajectory  (r\^  (t)  ,1^  (t) )  over  a 
planar  surface.  Aboard  this  object  is  a  dynamic  system  which  is  affected 
by  a  random  field  f •  He  would  like  to  consider  the  statistical  des¬ 
cription  of  the  state  x(t)  of  the  system  in  tens  of  the  specified  trajectory 
and  the  statistical  description  of  the  random  field.  Problems  of  this 
general  type  arise  in  applications  such  as  inertial  navigation  [1,2]  where 
f  represents  the  errors  in  our  knowledge  of  the  variations  in  gravity  and  x(t) 
consists  of  the  errors  in  an  inertial  navigation  system.  Since  the  inertial 
system's  accelerometers  measure  actual  acceleration  plus  gravity,  an  estimate 
of  gravity,  from  a  gravity  map  of  some  sort,  must  be  subtracted  from  the 
accelerosmter  outputs.  Thus  map  errors  directly  drive  the  dynamics  of  the 
navigation  system. 

In  the  next  section  we  develop  equations  for  the  evolution  of  the 
covariance  of  x(t)  for  a  particular  class  of  random  fields.  For  the  special 
case  of  a  straight  line  trajectory  that  reverses  on  itself,  we  develop  in 
Section  III  a  novel  Markovian  representation  for  the  process  x(t),  and  in 
Section  XV  we  use  this  representation,  together  with  recent  results  on  the 
real-time  updating  of  smoothed  estimates,  to  solve  an  estimation  problem. 


IX.  Covariance  Analysis  for  Motion  Through  a  Two-Dimensional 
Random  Field 


Let  f(n  ,rv)  be  a  two-dimensional  stationary  Gaussian  random  field 

1  2  ...  L-CJ’  \ 


which  for 


for  this  field  is 
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simplicity  we  assume  to  be  seto  mean.  The  correlation  matrix 
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E(f (t,s)f ' (0,0) ]  m  R(t,s)  (2.1) 

It  la  easily  seen  from  (2.1)  that 

R(t,a)  -  R* (-t,-s)  (2.2) 

Let  (H^(t) ,0^ (t) )  b«  a  specified  trajactory  through  tha  plana  and 
conaidar  a  dynamic  systam  driven  by  tha  fiald  along  tha  trajactory 

x(t)  -  Ax (t)  +  f  (n^t)  ,n2  (t) )  +  w(t)  (2.3) 

where  w(t)  la  a  zero  naan  white  Gaussian  procasa  with 

E[w(t)W  (T)l  «  Q<S(t-x)  (2.4) 

We  aaauaa  that  tha  Initial  condition  x(0)  ia  zaro  maan  and  Gauasian  and 
that  x(0),  w ,  and  f  ara  mutually  indapandant.  Wa  would  Lika  to  determine 
tha  avolution  of 

P(t)  -  E[x(t) x' (t)l  (2.5) 

wa  will  put  furthar  rastrictions  on  tha  fiald  f  that,  as  wa  will  see, 
laad  to  P(t)  baing  spaciflad  by  a  finita  sat  of  matrix  diffarantial 
aquations.  Spacifically  wa  will  assums  that  tha  covariance  R  is  saparabla 

R(t,T)  -  R1(t)R2(a)  (2.6) 

whara  wa  assume  that  R^  and  R2  ara  square  and  that 

R^t)  -  R^(-t)  ,  R^s)  -  R^(-s) 


(2.7) 


L,2 
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ttd 

Pit 

Ri(t)  "  Hi*  Gi  '  1-1 

(2.1),  (2.6)  and  (2.7)  wa  can  also  daduca  that  R^  and  coauta  Cor 
any  valuaa  of  thair  arqunanta .  this  nodal  la  tha  continuoua-tiaa  analog 
of  tha  nodal  ananinad  by  Attaai  [3].  Spacifically,  tha  2-D  apact ram  of  f 
ia  aaparabla  and  rational. 

Aa  a  firat  atap  in  obtaining  tha  daairad  aqua ti ana  for  P(t),  dafina 

Q(t,a)  -  R(n1(t)-nl(a) ,  n2(t)-n3(«)i  (2.9) 

Than,  writing 

t 

*(t)  -  aAt*(0)  ♦  J  aA(t"T)  [f(n1(T),na(T))+w(T))dT  (2.10) 

0 

wa  can  obtain  and  a^raaaion  for  P(t)  from  (2.5).  Diffarantiating  wa 
obtain  tha  baaic  aquationa 


P(t)  -  AP(t)  ♦  P (t) A'  ♦  L(t)  +  L‘ (t) 


(2.11) 


L(t) 


Q(t,T)aA' (fe"T)dT 


(2.12) 


Tha  problan  than  baconaa  ona  of  dataxaining  a  aat  of  diffarantial  aquationa 
for  L(t) .  This  calculation  da pan da  upon  tha  natura  of  tha  trajectory . 

Thara  ara  aavaral  caaaa  to  ba  axaninad.  For  simplicity,  wa  will  aaaun 
throughout  that  0^(0)  ■  n2(0)-0. 


Cui  li  This  is  ths  sinplsst  csss  in  which  vs  don't  changa  quadrants 
in  which  vs' rs  hsading .  That  is,  if  ws  chooss  ths  northsast  as  ths 
diraction  of  notion ,  ws  haws  ths  situation  dspictsd  in  Fiqurs  2.1a  whsrs 


nl<t>  "  nl<a)  -  0 


n2(t)  -  nx(s)  >  o 


(2.13) 


In  this  casa,  using  (2. 6) -(2. 9)  ws  find  that  (2.12)  can  ba  writtan  as 


f  v 


p, (n. (t)-n. (o))  f  (a,(t)-n, (a))  t  . 
11  l  „  „  _  2  2  2  _  a  (t-». 


G1H2* 


^  H^tt) 


(2.14) 


Diffsrsntiating  B^(t),  ws  obtain 


Bl(t)  -  r^tOFjB^t)  ♦ 


♦  n2(t) 


/fc  f  (n.(t)-n.  (o))  f  (n,  (t)-n„  (o)> 

s  x  1  <iiH2T2*  A(t-0)d0 

0  ... 


(2.15) 


Nota  ths  F^  factor  in  ths  aiddls  of  ths  last  tarn  of  (2.15).  This  lsads  to 
ths  following.  Da fins 

^t  f,  (n.  (t)-n,  (a))  .  ,  F,(n„(t)-n,(o))  .... 

IjW  •f  s  1  1  1  e^jFj  ls  2  2  2  «32sa  (t‘0,d0 


(2.16) 
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We  know  that  there  is  an  integer  r  and  coefficients  pQ, . . .  ,p^_^ 


9UU1  Ul«b 


r-1 


j 


f2  *  Z  »/2 

*  j-0  3  * 


(2.17) 


Than 


L(t)  -  (t) 

Bj(t;)  -  (tjFjBj  (t)  +  Bj  (t)A'  +  n2( t)Bj+1(t)  +  g1h2f^“1g2 

icj<r-l  (2.19) 


r 

B  (t)  -  f)  (t)F,B  (t)  +  B  (t)A'  +  fl  (t)  l  p  B 

r  i  j 


(t) 


* wrs 


(2.20) 


Bj (0)*0  ,  j-l,...,r  (2.21) 

Note  that  we  can  obtain  analogous  aquations  with  the  roles  of  F^  and 
reversed  if  we  use  the  comnuta tivity  of  R^(t)  and  R., (t) .  Thus  in  this 
case  we  obtain  a  finite  set  of  linear  matrix  differential  equations  for 
L  and  therefore  for  P.  Note  that  if  the  trajectory  is  a  straight  line 
—  i.e.  n^(t)«a,  n2 (t)-0  —  then  these  equations  are  time'- invariant  and 
are  equivalent  to  one  higher-dimensional  Lyapunov  equation  for  x  and  for  the 
state  of  a  shaping  filter  for  f  along  this  line. 


Case  2:  In  this  case,  illustrated  in  Figure  2.1b,  we  have  a  change  of 
quadrants  from  northeast  to  southeast.  Clearly  the  following  analysis  also 
holds  for  any  turn  from  one  quadrant  into  an  adjacent  one.  Mathematically, 
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nL(t)  -  n1<s) 

1  0 

t>s 

Hjtt)  -  n2(s) 

>  0 

*1*0 (t)  (2.22) 

Ojtt)  -  n2(s) 

i  0 

•*o (t) 

(t)  is  defined  in  the  figure. 

Here  t^  is  the  time  at  which  our 

tura  takas  us  into  snothsr  quadrant  in  direction ,  and  t_,  is  ths  tins  at 
which  (t) -0 . 

For  t£t^,  ths  analysis  of  this  cast  is  identical  to  that  for  Cass  1. 
Thus,  consider  t^<t<t2  and  1st  us  break  up  ths  integral  expression  for 
L(t)  s 


/•o{t>  f,  (n,  (tj-n.  (si)  r,(n,(t)-n,(s)) 

1  1  1  GlH2e  2  2  2  G2eA  #,ds 


♦ 1  v 


*  f,  (n.  (t)-n,  (s))  f' (n,  (s)-n.(t))  .....  . 

...11  1  - '22  2  „•  A(t-s) 


G1G2* 


Hje  'ds 


sQ(t) 


(2.23) 


where  we  have  used  the  fact  that  R2 (n^ (tj-n^ (s) )  ■  (a, (s)-n2 (t) ) .  in 

differentiating  (2.23)  we  will  need  to  calculate  sQ (t) .  This  can  be  done 
as  follows.  By  definition 


n2(sQ(t))  -  n2(t) 


(2.24) 


Therefore 


r\2(tQ(t))sQ(t)  »  n2(t) 


(2.25) 
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or 

n2  (t) 

«Q(t)  -  -  (2.26) 

Vo(t)) 

Not*  that  if  n(Sj(t))«0,  aa  it  is  in  Figure  2.2,  w*  will  havw  to  *valuata 
higher  derivatives.  This  causes  no  conceptual  difficulty  but  it  simply 
complicates  the  development.  Therefor*  we  will  assume  for  simplicity  that 
there  are  no  inflection  points  in  the  trajectory  over  the  interval  [0,t^) . 
Let 


s  (t) 


Bj(t)  -  y°  /i(Tii(t)-^(s))W2j-i/2^(t,^2(9))vA-(t-.] 


j-l,...,r 


(2.27) 


Note  that  if  w*  define  sQ(t)«t  for  t<t^,  then  B^  is  precisely  the  quantity 
in  equation  (2.16)  and  thus  the  initial  condition  at  time  t^  for  in 
(2.27)  is  Bj (t^)  calculated  from  (2.19)-(2.21)  .  If  we  now  differentiate 
(2.27)  and  use  (2.17)  w*  find 


Bj  (t)  -  n^tlF^lt)  ♦  Bj  (t)A'  +  fl2(t)Bj+i(t) 


^(t) 


F1(n1tt)-n1(sQ(t))]  A((t-s  (t)> 

G1H2F2  G* 


j-1, . . . ,r-l 


(2.28) 


r 

Br(t)  -  n^C t)FlBr<t)  ♦  Br(t)A’  ♦  n2(t)  l  p^B^ 


(t) 


n2(t)  F1[n,(t)-n1(s0(t))i  a* (t-s  (t)) 

G1H2F2  G26 


+  V90(t))  * 


ds 


(2.29) 


Note  that 


Bj  (t2)-0 
sine*  sQ(t2)-0. 


-11- 


(2.30) 


Now  let 


F  (n.  (t)-n  (3))  ‘  .  .  Fl(n  (s)-n  (t) )  ...... 

6  GlG2(F2)j  6  H’eA  (t"8) 


s0(t) 


ds 
(2.31) 


Then 


Cj (t)  -  n1(t)F1cj(t)  +  (t)A*  -  n2(t)cj+1(t) 


^  ,  n,(t)  f,  (n,  (t)-n,  (s_(t)))  .  . 

+  GlG2(F2)  H2  ■  n2  (s0 (t) )  6  G1G2(F2}  H2e 


A'  (t-S( 
(2.32) 

'  l-r,'  \  r~^Tr' 


”r(t)  "  Vt)Ficr(t>  +  cr(t)A'  -  n2(t)  l  p^c.  (t)  +  GjG^(f2)  h2 

n  (t)  f  (n.  (t)-n. (s  (t) ) )  A'(t-s.(t)) 

* _  _  1  X  XU  _  _  i  ,_i.r— X  u 

■  n2(sQ(t))  *  G1G2  F2  H2  (2.33) 


cyt^-o  , 


j«l, . . .r 


(2.34) 


Then 


L(t)  -  H^B^t)  +  Cx(t)  ] 


*1^2 


(2.35) 


(t) ) 


Note  that  in  the  case  of  a  piecewise  linear  trajectory,  such  as 
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\(t)  -  ! 

i“» 

i 

t<fcl 

i 

a-> 

t>t. 

1 

1 

t<t 

Vt>-  j 

1  1 

t>t, 

•0<t)  - 


K 


(2.36) 


(2.37) 


(2.38) 


vhich  la  nagativa  hara  ainca  (^>0,  B2<0. 

Wa  now  naad  only  piaca  togathar  tha  aituation  for  t>t2>  In  thia  caaa 


,,  ,  /.fc  F  (n.(t)-n  (a))  r’(n2(a)-n,(t)) 

L<t)  -  f  H^a  1  G^a2  2  2  hJ.* (2. 


39) 


Thus  in  thia  ragion 


L(t)  -  HjC^t) 


^2 


whara  Cj^tj)  ia  obtainad  from  (2.32)-(2.34) ,  and,  for  t>t, 

C^t)  -  Al<t)rlcJ(t)  ♦  c(<t)v  -  tHj(t)cj+1(t) 


♦  <31<32(F2)3'1h2 


(2.40) 


J»l, . . .r-1 


(2.41) 


13' 


.  t 

-r  {t)  "  Vt)FlCr(t)  +  cr(t)A’  -  ^2 (t)  l  Pj-lCj 


(t) 


* 


(2.42) 


Thus  in  Cass  2,  ovsr  the  time  interval  [0,t^]  we  have  one  set  of 
equations  to  calculate  L;  over  the  interval  [t^,t2] '  whi^®  th®  n2  coordinate 
of  the  trajectory  retraces  its  path  over  [0,t^],  we  have  two  sets  of  equations; 
and  for  t>t^  we  are  back  to  a  single  equation  which  is  essentially  the 
equation  obtained  in  Case  1,  except  here  we  have  a  southeasterly  trajectory 
as  opposed  to  the  northeasterly  trajectory  of  Figure  2.1a.  It  should  be 
clear  that  we  can  do  this  for  arbitrary  trajectories.  Only  during  the 
"transient"  of  a  turn  do  we  pick  up  additional  equations.  In  the  Appendix 
we  describe  one  somewhat  more  complex  case  in  which  both  and  n2  coordinates 
simultaneously  retrace  previous  values  (this  does  not  mean  a  trajectory  that 
retraces  itself  —  see  Figure  A.l) .  In  that  case,  there  are  two  additional 
sets  of  equations.  From  the  cases  considered  in  this  section  and  in  the 
Appendix  it  is  not  difficult  to  see  that  at  any  time  we  must  include  m  additional 
sets  of  equations,  where  m  is  the  total  number  of  previous  times  in  the 
trajectory  that  either  the  r|^  or  n2  coordinate  of  the  trajectory  equals  the 
corresponding  coordinate  at  the  present  time.  If  A  is  a  stable  matrix,  then 
the  effect  on  L(t)  (and  hence  P(t))  of  a  trajectory  turn  far  in  the  past  becomes 
insignificant.  This  can  be  seen  in  (2.28)  where  the  driving  term  goes  to  zero 
exponentially  as  t-sQ(t)  -*■«  (the  matrices  F^  and  F2  are  stable  since  f  is  a 
stationary  process  with  finite  covariance) .  Thus  in  practice  we  need  only 
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Jc««p  track  of  turn*  within -a  cartain  number  of  time  constant*  of  A  and 
correlation  distance*  of  the  field  (inverses  of  the  magnitudes  of  the 
eigenvalues  of  and  ?2) . 

Hi.  Markov- Type  Models  for  Over-and-Back  Trajectories 

In  the  preceding  section  and  in  the  appendix  we  performed  some 
relatively  straightforward  calculations  to  obtain  sets  of  differential 
equations  for  the  propagation  of  the  covariance  of  the  state  of  a  dynamic 
system  moving  through  a  random  field.  The  primary  contribution  of  that 
analysis  is  to  provide  some  understanding  of  how  the  geometry  of  the  tra¬ 
jectory  affects  the  state  covariance.  While  this  is  of  some  use,  there  is 
sti11  •  great  deal  left  to  understand  about  the  fundamental  way  in  which 
the  uncertainty  in  the  field  affects  the  statistics  of  the  process  x(t). 

In  this  section  we  will  develop  a  Markov- like  description  for  the  special 
case  of  over- and- back  trajectories.  This  not  only  provides  us  with  further 
insight  into  the  evolution  of  x(t)  but  it  also  forms  a  basis  for  solvir~ 
estimation  problems  of  this  type,  a  topic  which  is  considered  in  the  next 
section . 

The  case  that  we  will  examine  in  this  and  in  the  next  section  involves 
a  trajectory  consisting  of  a  straight  line  path  followed  by  a  reversal  of 
direction  and  a  return  trajectory  over  the  same  path.  W*  also  will  assune 
a  constant  velocity  (normalized  to  1)  over  both  segments  of  the  path,  but 

assumptions  is  made  only  for  clarity  in  our  exposition  as  is  our 
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assumption  that  tha  dynamic  syatam  is  time- invariant. 

It  should  ba  claar  from  our  analysis  as  to  how  our  rasults  can  ba  modifiad 
to  account  for  ncnuniform  velocity  and  time- varying  systems. 

Consider  tha  nodal 


i(t)  -  Ax(t)  +  u(t) 

♦  w^(t) , 

0<t<2T 

(3.1) 

where  w^  is  a  white  noise  process  with 

E(w^ (t)w^  (t)  ]  -  S^ft-T) 

(3.2) 

and  where 

(  f(t) 

0<t<T 

u(t)  »  \ 

(3.3) 

[  f  <2T-t) 

T<t<2T 

Here  f  is  a  one-dimensional 

process  (representing  the  field  along 

the 

track,  and  we  assume  the  f  can  be  modeled 

as  the  output  of  a  finite- 

dimensional  shaping  filter 

* 

5<t)  -  F5(t)  +  w2<t) 

0<t<T 

(3.4) 

f(t)  -  H?(t) 

(3.5) 

where  v>2  is  white  noise,  with 

E tw2 (t)w£ (t) ]  -  S25(t-T)  (3.6) 

We  assume  that  all  of  the  processes  above  are  sero  mean  and  Gaussian  and 


that  x(0) ,  w^,  5(0),  and  w2  are  mutually  independent. 
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Por  0<t<T  wa  have  tha  sum  situation  as  In  Casa  1  consldsrad  in  tha 
preceding  s action.  Over  this  tins  interval,  while  we  ara  going  forward, 
tha  joint  procass 


«2(t)  - 


C(t) 

x(t) 


(3.7) 


is  Markovian ,  with  tha  following  stata  aquation 


*2(t)  - 


P  0 

H  A 


x2(t)  + 


w2  (t) 


LWl(t> 


(3.8) 


0<t<T  (forward) 


Tha  manning  of  tha  notation  (forward)  in  (3.8)  will  bacoma  claar  shortly. 
Thus  tha  covarianca  L (t)  of  *2(t)  can  ba  obtainad  from  tha  diffarantial 
aquation 


"p  0 

'  F*  H'" 

— 1 

CO 

J 

o 

— 1 

L<t>  - 

L<t)  +  L(t) 

+ 

«>  1 

2 

H  A 

2  2 

1 

o 

> 

1 _ 

1 — 

o 

CO 

(3.9) 


A*  w*  saw  in  Casa  1,  wa  naaded  ona  additional  sat  of  aquations  in  ordar  to 
calculata  tha  covarianca  of  x  (t(L  Hara  wa  saa  that  that  sat  of  aquations  as- 
s antially  comas  about  by  augmanting  tha  stata  x(t)  with  a  shaping  filtar 
modal  for  tha  fiald  in  ordar  to  obtain  a  procass  that  is  Markovian.  Ones 
that  is  dona,  as  in  (3.8),  wa  can  usa  standard  results  to  writs  down  tha 
covarianca  aquation  (3.9). 


I  J 
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The  interesting  pare  of  this  analysis  occurs  ovar  the  time  interval 
T<t<2T,  since  here  we  are  reversing  over  the  sane  sample  path  of  £.  Again 
our  goal  is  to  augment  x(t)  with  something  in  order  to  obtain  a  Markov  model 
over  this  time  interval.  In  order  to  do  this  we  clearly  must  consider  a 
model  for  £  that  runs  in  reverse.  Using  the  results  in  [4]  we  can  write  a 
reverse  time  modal  for  the  augmented  process  x^ (t) : 


- 

•  w 

_ 

P  0 

S2  ° 

-1 

w2  (t) 

+ 

A 

l  <*> 

x  (t)  - 

2 

2 

w,  (t) 

‘ 

H  A 

-°  Sl. 

L  1  J 

0£t£T  (backward)  (3.10) 

Here  (w^ (t) ,  w^(t))  is  a  white  noise  process  backward  in  time  independent  of 

xj (T)  ■  (£* (T) ,x' (T) ) .*  The  processes  have  the  same  statistics  as  the 

w^.  Note  one  interesting  aspect  of  this  model.  If  we  examine  (3.8)  or 

(3.4)  we  see  that  £  is  a  Markov  process  by  itself  forward  in  time  —  i.e. 

it  is  decoupled  from  x(t).  However  this  is  not  true  in  the  reverse-time 

P-1 

modal  (3.10)  since  is  not  block  diagonal.  The  reason  is  that  going 

forward  in  time  the  process  £  drives  the  process  x.  Then,  since  the  reverse 
process  is  a  Markovian  representation  of  x2  given  its  future,  we  should 
expect  to  see  coupling,  since  the  present  value  of  £  is  certainly  not  independent 
of  the  future  of  x. 

*  The  use  of  white  noise  here  makes  our  derivation  somewhat  informal.  However 
it  is  conceptually  correct.  To  be  precise  we  should  replace  w^  by  dfi^  where 

the  8^  are  Brownian  motion  process e  forward  in  time .  Than  w^  is  replaced  by 

d6^  which  is  a  Brownian  motion  backward  in  time ,  independent  of  xj (T)  « 

xj (T)  -  (£' (T) ,x^(T) ) .  This  process  is  obtained  by  subtracting  from  d6^  that 

part  which  is  predictable  given  the  future  of  .  See  [4]  for  details. 
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Zf  we  new  let 


«3(t)  -  Xj  C2T-t) ,  nA(t)  -  w^T-t) 
we  obtain  a  nodal  forward  in  tlma  ovar  tha  Interval  T<t<2Ts 


(3.11) 


j 

m  m 

P  0 

I 

(A 

KJ 

o 

_ 1 

-1  ) 

n2  (t) 

£  (t)(  X  (t)  + 

( 

H  A 

1 

O 

M 

*-• 

1 _ 

2  )  3 

.  ni(t)  - 

T<t<2T  (forward) 


(3.12) 


Tha  initial  condition  for  this  process  is  *3 (T)  -  x2 (T) ,  with  covariance 

l  (t). 

2 

Consider  now  tha  following  augmented  process  over  the  tine  interval 
T<t<2Ts 


*4 ft)  - 


x3(t)  ‘ 

’  5(2T-t)  ' 

- 

x(2T-t) 

x  (t ) 

.  x(t ) 

(3.13) 


Then,  using  (3.1),  (3.3),  (3.5),  and  (3.11)  we  obtain  a  Markovian  repre¬ 
sentation  for  the  behavior  of  this  augmented  state 

.-1 


.  r-io  ; 

L  (0,  H) 


(t) 


AJ 


rn2(t) 

x.  (t)  +1  nl(t) 
4  1  wx(t) 


(3.14) 


-19- 


where 


n2  (t) 

- 1 

u 

o 

o 

-i 

E 

nL(t) 

V  (t) 

Jnj(T) ,n^(t) ,u^(T))j  - 

0  0 

0  0  S, 

i 

L 

-  1  J 

«(t-T) 


Basically  (3.14)  describes  a  method,  starting  at  t»T,  for  a iaul tanaoua ly 
generating  tha  futura  (t>T)  of  x(t)  and  its  past  (t<T) .  In  this  fashion 
va  can  taka  into  account  tha  fact  that  tha  trajactory  has  ravarsad  its 
diraction. 

Ms  can  usa  (3.14)  as  tha  basis  for  detexaining  tha  covariance  for  x(t) 
Spacif ically ,  dafina  N(t)  as 


M(t) 


•![:  :]*[:•  :J  c 


(t)( 


(0,  H) 


Than,  latting  £  (t)  danota  the  covariance  of  x.  (t),  we  obtain 
4  4 


l  (t)  -  N(t)J  (t)  +  l  (t)M'(t)  + 

4  4  4 


s2  0  0 

0  0 

0  0  s 


1 J 


(3.15) 


(3.16) 


To  obtain  tha  initial  condition  for  this  equation,  note  that 


5(T) 
x  (T) 
x(T) 

■ 


X^(T)  • 


(3.17) 
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Thus,  if  w«  writs 


/Mu  (VT,)u\ 

\(J2(T1)u  fe»W 


.then 


(VT,)u  (V")u\ 

M**  (V*^) 

(l2<«)22  (I2<T>)2/ 


(3.18) 


(3.19) 


Thus  ws  sss  that,  as  in  Case  XI,  a  reversal  of  motion  leads  to  an  ad¬ 
ditional  equation.  Also,  one  can  regard  the  over-and-back  example  as  a 
degenerate  form  of  Cases  b  and  c  which  are  examined  in  the  Appendix 
(referring  to  the  notation  in  the  Appendix,  in  the  over-and-back  case 
s^(t)  •  s2 (t)  for  all  t  and  t^  “  t^).  Thus  the  straightforward  analysis 
of  the  Appendix  will  lead  to  equivalent  equations  in  this  case. 

Note  that  based  on  the  understanding  gained  in  this  and  in  the  preceding 
section,  we  can  see  what  will  happen  for  more  general  over-and-back  tra¬ 
jectories.  For  example,  as  illustrated  in  Figure  4.1a,  consider  the  case  in 
which  we  continue  the  process  for  t>2T  without  any  further  change  in  course. 
Xt  is  not  difficulty  to  show  that  for  t>2T  we  cam  once  again  obtain  a 
Markovian  representation  for  the  joint  process 


where  the  initial  covariance  £ (2T)  for  this  process  at  timm  jt  is  obtained 

* 

fro*  the  solution  to  (3.16) 

(V”’)u  (V”')u 

In  this  case  the  tine  period  [T,2T]  represents  a  transient  due  to  the  turn , 
whose  effective  will  becone  negligible  if  A  is  seable.  In  fact  in  this 
case  as  t-*«,  x  will  achieve  the  same  steady-state  covariance  in  this  situation 
as  one  would  fro*  a  trajectory  that  moves  to  the  left  for  t>0  without  any 
turns.  Similarly,  if  we  consider  a  second  course  reversal  as  in  Figure  4.1b, 
we  must  obtain  a  reverse  time  model  for  x.,  reverse  time  once  again  to 
obtain  an  equation  for  x^(t)  «  x^(4T-t)  and  augment  this  with  x(t)  to  obtain 
a  Markovian  model  over  the  time  period  2T^t<3T.  Thus  in  this  case  we  obtain 

r  t . 

another  additional  equation  for  the  covariance  evolution. 

IV.  Over- and- Back  Estimation 

In  this  section  we  consider  the  problem  of  estimating  the  process 
described  by  (3.1)-(3.6)  given  measurements .  Specifically,  suppose  we  assume 
that  the  random  field  has  been  mapped  by  a  previous  survey 

y^(t)  -  Cjt(t)  ♦  v^(t) , 


0<t<T 


(4.1) 
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where  E  [v^  (t)  v^  (T)  ]  ■  R^6(t-T),  so  that  wa  hava  the  smoothed  estimates 

£#(t)  -  E[£(t)|  y^T),  0<t<T]  (4.2) 

Consider  now  a  set  of  seal-time  measurements 

*(t)  -  C2*(t)  +v2(t),  0£t£2T  (4.3) 

where  E(v2  (t)v2  (t)  1  -  R2<$(t-T)  and  v^  and  v2  are  independent.  We  wish  to 
consider  the  problem  of  using  the  previously  mapped  information  (4.2), 
together  with  the  new  data  (4.3)  to  estimate  x(t). 

As  in  the  preceding  section,  this  problem  is  best  analyzed  by  considering 
the  two  intervals  [O.T]  and  [T,2T]  separately.  Thus,  let 

y2(t)  -  z(t),  te(0,T],  y3(t)  -  z(t),  te[T,2T]  (4.4) 

The  problem  over  the  first  time  interval  is  a  real-time  smoothing  problem, 
that  is,  we  have  smoothed  estimates  for  part  of  the  state  (here  the  £  part 
of  *2  as  given  in  (3.7)  and  (3.8))  from  previously  taken  data,  and  wish  to 
incorporate  these  estimates  into  an  overall  state  estimate  given  new  real- 
tiae  data.  A  real-time  smoothing  problem  of  this  type  was  solved  in  [5] . 

In  order  to  apply  that  solution  here,  define 


X2(t|T,t)  -  E(x2 (t)|  y^ (t) ,  0<T<T;  y2 (t) ,  0<T<t] 

Then 

x2(t|T,t)  -  P2(t) (P^1(t)qf (t)  +  P^1(t)qb(t)] 


(4.5) 


(4.6) 
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where 


^(t)  - 

qf (0)-0 

-%(t)  " 
q.(T)-0 


'o 

A 

’  0 

Ff(t)qf(t)  + 

C#(t)  +  Pf(t) 

-1 

H 

o 

K» 

J* 

l— 

Fb(t)qb(t)  +  (Pb(t)  l  (t)-I) 


(4.7) 

(4.8) 

(4.9) 

(4.10) 


when  ^(t)  is  the  unconditional  covariance  of  (t) ,  as  calculated  from 
(3.9).  The  remaining  quantities  in  (4.5)-  (4.10)  are  deterministic  and 
are  determined  from  the  equations 


F  0 

F'  H* 

S,  o 

• 

P,(t)  +  P.(t) 

+ 

H  A 

E  E 

0  A' 

°  Sl- 

-  Pf(t) 


0 


C.R~^C_ 
2  2  2 


Pf  (t) 


(4.11) 


Pf  (0) 


■Pb(t) 


1,(0) 


F  0 


H  A 


S2  ° 


S1J 


,-l 


1 


l2  <*>  |V*> 


(4.12) 


Pb(t) 


F'  H' 
0  A* 


--1 

+  I  (t) 
2 


S2  ° 
0  S 


). 


-  V*’ 


ci\lci 


0 

0 


1J 


Pb(t) 


S2  ° 
0  s 


1 J 


(4.13) 
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P.  (T)  -  I  (T) 
b  2 

P^tfe)  -  P'l(t)  +  P*l(t)  -  l  1(t) 

2  r  b 


F  (t)  -  -  P  (t) 

r  H  k 


C  *  R~  0 

11  1 


CiR2lc2 


Fb(t) 


F  °  S2  °  ! 


H  A  0  S 


L  (t)  -  p„(t) 


c '  R~ 'c  0 

11  1 


0  0 


(4.14) 


(4.15) 


(4.16) 


(4.17) 


Let  us  briefly  interpret  these  equations.  For  a  detailed  discussion, 
we  refer  the  reader  to  [5].  Here  Pf(t)  is  the  estimation  error  covariance 
for  x2(t)  given  y^(T)  and  y2(t),  T£t  --  i.e.  the  causal  estimate  of  x2(t) 
using  only  causal  information  from  the  previously  collected  data  y^  and  the 
new  data  y  .  Similarly,  P.  (t)  is  the  estimation  error  for  a  reverse-time 
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filter  which  estimates  x2(t)  given  only  y^T),  t<T<T.  The  set  of  information 
used  in  these  two  estimates  comprise  all  the  information  used  in  computing 
x,(t|T,t).  The  forward  filter  (resulting  in  an  estimate  with  covariance 
Pf(t))  is  the  usual  Kalman  filter  and  has  the  form 


xf(t)  -  Ff(t)xf(t)  +  Pf(t>r  c’r" 


0  YL(t) 
C2R21  y2(t) 


(4.10) 


Similarly  the  backward  Kalman  filter,  with  estimation  error  covariance 


P.  (t) ,  has  the  form 

D 
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Pb(t>*b(t>  +  Pb(t) 


(4.19) 


lb*  estimation  error  covariance  for  x^(t|T,t)  is  P^  (t)  ,  and  the 
expreasion  for  it  in  (4.15)  ie  taken  from  Mall  [6].  In  fact,  using  the 
smoothing  equations  in  [6]  and  the  result  of  (5),  x2 (t|t,T)  can  be  written 


Xj  (t  |T,t)  -  p2(t)  tPjl(t)xf(t)  ♦  P^ttJx^t)]  (4.20) 

The  fact  that  £  ^ (t)  is  subtracted  on  the  right-hand  side  of  (4.15)  reflects 
the  fact  that  the  estimates  xf(t)  and  ^(t)  of  the  state  are  correlated,  as 
they  both  use  the  a  priori  information  which  has  an  uncertainty  specified  by 
the  unconditional  statistics  of  x^ (t) . 

Finally ,  (4.6)-(4.10)  are  obtained  from  (4.18)-(4.19)  by  some  manipula¬ 
tions  aimed  at  replacing  y^(t)  in  (4.18),  (4.19)  by  the  previously  determined 
map  £s<t).  The  details  of  these  calculations  in  a  somewhat  more  general 
context  can  be  found  in  (5).  Note  that  only  q^(t)  is  driven  by  the  new,  real- 
time  data  y„(t),  while  q.  (t)  is  a  functional  only  of  the  smoothed  map  E  (t)  and 
in  principal  can  be  precomputed.  In  practice,  what  this  means  is  that  once 
the  one-dimensional  trajectory  has  been  charted,  we  can  integrate  (4.9) 
backward  along  this  trajectory,  store  the  result,  and  combine  it  with  the 
stored  map  E^(t)  and  q^(t)  through  (4.6)  in  order  to  determine  x,(t|T,t). 

Consider  now  the  estimation  problem  over  the  time  interval  (T,2T). 

This  is  again  a  real-time  smoothing  problem,  thanks  to  our  augmented 
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Markovian  modal  (3.13),  (3.14).  Specifically  using  tha  data  {y^t),  y2 (T) } 
wa  hava  conputad 

x3#(t)  -  E(x3(t) |yx (x) ,y2(x) ,  0<T<Tl 

-  B[x2(2T-t)  y^T),  y2  (t)  ,  0<t<T]  -  x2>(2T-t)  (4.21) 

Tha  quantity  that  wa  wish  to  compute  is 

x4(t|T,T,t)  -  E[x4(t) |yi(T),  0<T<T,  y2(T),  0<t<T,  y3(t),  T<t<t] 

-  E[x4(t) |yx (T) ,  0<T<T,  2(T),  0 <T <t ]  (4.22) 

Using  tha  raal-tlaa  smoothing  formulas  of  [5] ,  wa  can  axprass 

x4(t|T,T,t)  in  tanas  of  x^  and  y3  for  te[T,2T) .  However,  thara  is  ona 

complication  caused  by  tha  fact  that  £  (T)  is  not  invertible  (see  (3.19)). 

4 

This  is  due  to  tha  fact  that  x(2T-t)-x(t)  for  t-T.  In  order  to  make  tha 
necessary  modifications,  it  is  convenient  to  change  basis.  Let 
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Than,  fro*  0.14),  (4.23),  and  (4.24) 


Furthar,  lat  us  assum  that  x(t)  as  dafinad  in  (3.1)  is  controllable  froai 
tha  noisa  w^(t) ,  i.a.  that  (A,S^)  is  a  completely  controllable  pair.  In 
this  case  it  is  not  difficult  to  sea  that  for  any  t>T,  V(t)>0.  than, 
da fining 
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P  (t  |  T,T,  t)  - 


ZOO 
OIO 
0  -I  I 


X4 (t|T,T,t) 


(4.28) 


w«  can  adapt  the  reaults  in  [5]  to  obtain 


p(t|T,T,t)  -  E(t)^(t)rf(t)  +  M^(t)rb(t)J  +£*3*(t)j 


(4.29) 


rf(t)  "  6f(t,rf(t)  + 


(-T2i(t)rH+A-T22(t)) 


x3#(t)+Mf (t) 


c'r’1 

C2R2 


y3(t) 


(4.30) 


rf (T)-0 


(4.31) 


-  Gb(t)rb(t)  +  (Hb(t)v',’(t)-i) 


("T21(t),H+A'T22(t)) 


*3,(t) 


(4.32) 


rb(2T)-0 


(4.33) 


where  M^,  E,  G^,  and  G^  are  given  by 


M  (t)  -  D(t)M. (t)  +  M. (t)D' (t)  ♦ 


s2  0  0 

0  8X 

0  -*i  2Si 


CWC1 


0  iCjR^C. 


L  0 


C2Rl\  C& 


.-’I, 


Mf(t) 


'2  -* 


(4.34) 
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Mf  (T) 


(4.35) 


VT(t)l  ^(t) 


-^(t)  Id*  (t)  +  vT(t) 


S2  .  o  0 

0  S1  -S1 

.°  -S1  2S1-I 


-v° 


’ci\lcl 


0 

L  0 


0 


0 

0  J 


S,  o  0 

2 

o  sx  -sx 

-c  -S1 25 1  J 
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^  (2T)  -  V4(2T) 


(4.37) 


E+(t)  -  M^(t)  + 

*£<*> 

-  v+(t) 

m 

Gf(t)  -  D(t)  -  Mf(t) 

ci\\ 

0 

0 

2C'r"1C. 

2  2  2 

0 

C^l°2 

G^tt)  -  -D (t)  - 

s. 

2 

0  0 

V+(t)-Mb(t) 

0 

sl^l 

0 

*1*1- 

(4.38) 
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C&S 


(4.39) 

0 

0 


0 


0 
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The  various  quantities  hare  play  analogous  rolas  to  thosa  playad  by  tha 
corrasponding  quantltlas  in  (4.6)- (4.17) .  For  axampla,  E(t)  is  tha  arror 
covariance  assoclatad  with  tha  estimate  p(t |T,T,t) .  Tha  only  difference 
in  this  casa  is  tha  usa  of  pseudo- inverses  for  V,  and  E.  Thasa 

all  raprasant  covariancas  of  p,  and  thus  from  (4.23)  at  time  T  thay  ara  all 
singular.  However  bacausa  of  our  noisa-controllability  assumption,  all  of 
thasa  quantities  ara  positiva  dafinita  for  t>T.  Furthermore ,  it  is  not 
difficult  to  check  that  tha  estimate  p  at  time  T,  as  defined  by  (4.29)  does 
have  its  last  block-component  equal  to  zaro. 

Finally,  let  us  comment  on  tha  issue  of  computing  x3g (t)  or,  equivalently 
x2a (t) ,  as  defined  in  (4.21).  Recall  that  what  wa  calculated  over  tha 
interval  [0  ,T]  was  tha  real-time  smoothing  estimate  x2  (t  |T,t) ,  as  defined 
by  (4.5) ,  and  using  results  from  (5]  wa  displayed  an  algorithm  for  performing 
this  calculation  in  terms  of  tha  new  data  y2  and  the  previous  smoothed  estimate 

A 

trajectory  £  .  Tha  problem  wa.  wish  to  solve  now  is  tha  calculation  of  tha 
s 

smoothed  estimate  x2a(t)  using  y2  and  This  is  a  problem  in  tha  updating  of 

smoothed  estimates,  which  is  also  examined  in  [5].  Using  thosa  results,  we 


find  that 


*2s(t)  "  p2,(t)[pi1(t)<lf(t)  +  ^^(t)]  ♦ 


$a(t> 
L  o 


(4.41) 


where  P^,  and  q^  a re  as  before,  and 


-1 


-q^t)  -  Fr(t)qr(t)  +  [Pr(t)£  (t)-I] 


0 

A 

‘  0 

H 

5#(t)  +  Pr(t) 

_i  -l 

-  C2R2  . 

y2(t) 


(4.42) 


q^Tj-o 


(4.43) 
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and  Pf,  P2#,  and  as*  defined  by 


(4.47) 

Comparing  (4.41)- (4.47)  to  (4.6)- (4.17)  we  see  that  the  only  change  in 
computing  the  full  smoothed  estimate  (t)  versus  the  real-time  smoothed 
estimate  x^ (t|T,t)  is  the  incorporation  of  the  new  data  (y2)  into  the 
re verse- time  processor  (4.42),  (4.47)  and  into  the  corresponding  reverse¬ 
time  error  covariance  P  and  the  overall  smoothed  error  covariance  P 

r  2s 


V.  Conclusions 


In  this  paper  we  have  examined  the  effect  of  a  random  field  on  a 
linear  dynastic  system  moving  through  the  field.  He  have  developed  a 
methodology  for  calculating  the  covariance  of  the  state  of  the  dynamic 
system  along  any  trajectory.  The  evolution  of  this  covariance  is  clearly 
dependent  upon  the  nature  of  the  trajectory,  and  our  results  indicate 
explicitly  how  this  dependence  is  reflected  in  the  differential  equations 
that  must  be  solved  to  determine  the  covariance. 

In  the  case  of  one-dimensional  motion  we  have  gone  several  steps  farther 
in  our  understanding  and  analysis  of  over-and-back  trajectories.  Specifically, 
with  the  use  of  the  technique  for  constructing  backwards  Markovian  models 
we  have  developed  Markov  models  over  each  separate  undirectional  segment 
of  the  trajectory.  The  dimension  of  these  models  decreases  when  the 
trajectory  goes  beyond  the  region  covered  in  previous  segments  and  increases 
when  there  is  a  turn.  Using  this  model  and  results  on  real-time  smoothing 
we  then  were  able  to  solve  an  over-and-back  estimation  problem. 

Several  directions  for  further  work  suggest  themselves.  The  first  is 
the  detailed  investigation  of  the  estimation  problem  discussed  in  Section 
IV.  While  we  have  described  the  solution  to  this  problems  we  have  not 
exploited  its  structure  as  fully  as  is  possible,  eithexf  in  terms  of  obtaining 
efficient  on-line  solutions  or  of  gaining  insight.  For  example,  it  is  clear 
that  the  measurements  of  the  state  x  of  the  dynamic  system  provide  information 
about  the  field  £.  How  is  this  information  incorporated  in  the  solution 
of  the  real-time  smoothing  problem?  This  is  potentially  important  in 
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problems  in  which  we  wish  to  use  the  dynamic  system  to  estimate  the  random 
field.  Gravity  mapping  using  inertial  instruments  is  a  potential  application. 

A  more  significant  extension  of  our  work  is  the  development  of 
Markov— type  models  and  the  corresponding  estimation  algorithms  for  more 
general  2-0  trajectories.  This  will  involve  a  significant  extension  of  the 
notion  of  a  backwards  Markov  process. 

Finally,  an  important  generalization  of  the  problems  considered  in 
this  paper  are  to  systems  moving  along  trajectories  which  are  random 
themselves.  Specifically,  referring  to  our  general  model,  suppose  that 
(ri^  (t) ,  n2  <t) )  are  in  fact  components  of  x.  This  is  in  fact  a  more  realistic 
model  in  seme  applications.  While  our  results  do  not  address  this  problem, 
they  may  be  of  value  in  the  case  in  which  the  trajectory  is  only  slightly 
disturbed  from  same  nominal.  In  that  situation  our  analysis  might  form 
the  basis  for  a  perturbation  analysis  of  the  random  trajectory  problem. 
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APPENDIX 


W«  consider  one  final  case  of  the  trajectory-covariance  problem  of 
Section  II.  In  contrast  to  Case  2  in  which  the  turn  takes  us  into  a 
neighboring  quadrant,  we  now  consider  a  trajectory  that  has  a  sharp  angle 
and  takes  us  into  the  opposite  quadrant  —  i.e.  northeast  (NE)  to  south¬ 
west  (SW) .  Three  cases  of  this  will  be  considered,  and  these  are 
illustrated  in  Figure  A.l.  As  can  be  seen,  these  cases  represent  suc¬ 
cessively  sharper  turns.  In  each,  t^  is  the  time  when  the  direction  of 
the  direction  of  the  trajectory  changes  from  NE  to  SE,  and  is  the  time 
we  change  from  SE  to  SW.  Also  denotes  the  time  at  which  the  trajectory 
crosses  the  02  axis  in  the  southwesterly  direction.  This  corresponds  to  the 
time  at  which  the  coordinate  of  the  trajectory  has  evolved  from  0  to  its 

maximum  value  n, (t_)  and  has  decreased  back  to  zero.  The  time  t  is 
13  y 

defined  in  an  analogous  fashion . 

The  distinguishing  characteristics  of  these  three  trajectories  are 
as  follows:  in  trajectory  (a)  only  one  of  the  coordinates  at  a  time 
retraces  past  values  —  first  n2(t)  and  then  n^t);  in  trajectories  (b)  and 
(c)  both  coordinates  retrace  put  values  over  a  common  portion  of  the 
trajectory;  in  case  (b)  n2(t)  completes  its  retracing  before  (t) ,  while 
exactly  the  opposite  is  true  in  cue  (c) .  As  we  mentioned  in  Section  II, 
we  will  see  that  at  any  time  we  get  one  set  of  additional  equations  for 
each  component  that  is  retracing  put  values.  Thus  in  cue  (a),  we  will 
have  one  additional  set  of  equations  beginning  at  t^  and  ending  at  t^. 
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A  different  single  set  of  equations  begins  at  t  and  end  at  t  .  In  case 

3  x 

(b)  we  have  an  additional  set  of  equations  beginning  at  t^,  a  second  set 
is  added  at  tine  t^ ,  the  first  set  ends  at  t^ ,  leaving  us  with  the  second 
set  which  ends  at  t^.  in  case  (c)  the  only  differences  are  that  there  is 
a  change  in  the  equations  at  time  tQ,  and  also  the  second  set  of  equations 
ends  first  (t  <t  ).  This  verbal  description,  together  with  the  following 
analysis ,  should  make  clear  the  approach  that  can  be  taken  in  analyzing  a 
general  trajectory. 

Case  a:  Referring  to  Figure  2.1b  and  Figure  A. la,  we  see  that  the  present 

case  is  identical  to  Case  2  considered  in  Section  II  up  to  time  t3>  That 

is,  for  O^t^t^  we  use  the  equations  of  Case  1,  with  L(t)  -  H^B^t)  (see 

(2.18)-  (2.21) ) .  For  tjCtftj,  L(t)  -  H^B^t)  +  C^t)]  (equations  (2.28)- 

(2.35)),  while  for  t  <t<t  ,  L(t)  -  H.C,  (t)  ((2.40)-(2.42) ) .  For  t<t<t 

2 - 3  11  3 - x 


Hj^t)  - 

n1(s)>  o 

for 

0<3<31  (t) 

J 

ft 

1 

n1(s)£  o 

for 

9]_ 

(A.l) 

n,(t)  - 

n,(s)<  o 

for 

0<s<t 

2 

Coopering  (A.l) 

to  (2.22)  we  see  that 

this 

is  quite  similar  to  Case  2, 

except  here  we  i 

ire  turning  from  SE  to 

sw. 

Thus,  in  analogy  with 

Case  2 

we  can  write 
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f.  (n,  (t)-n,  («))  f'  (n,  (»>-n,  (t) )  ...... 

11  1  g.g' •  2  2  2  h>a  (t_,)d . 

12  2 


/f' (n.  (•)-n.  (t))  f' (n.(«)-n2(t)) 

V  "A*  V  * 


•  ^(t} 


(A.  2) 


I*(t)  -  (t)  G^D  (t> 


(A. 3) 


whsrs 


VV 

°j (t3> 


valuss  calculated  using  (2.40)-(2.42) 


(A. 4) 


c^t)  -  n^tiF^tt)  +  (t)A’  -  n2(t)cj+1(t) 

n,  (t)  .  f'  (n,  (s.  (t))-n,  (tn  a’ (t-s.(t)) 

G1G2(F2)  *  H2® 


j»l, . .  .r-1 


(A. 5) 


cr(t)  -  n 


L. 

Vt>FlCr(t)  *  cr(t)A*  *  na<t)  l  Pj.iCj 


n.  (t)  .  Fl(n,  (s  (t))-n-(t»  a*  <t-s  (t>) 

+  (*x  Ct) )  G1G2(F2)  •  H2 


(A. 6) 


Note  that  C.(t  )-0.  Also 
3  * 


Dj(t)  -  -^(tJP^tt)  +  Dj  (t)A'  -  H2(t)Dj+1(t)  + 


n.(t) 

^(tn  h1g2(f2 


A' (t-a,  (t) ) 


.  .  f  (n, (s  (t))-n, (t) )  “  "■  "1 

:G!(F:)jle2  21  2  v 


j»l, ...,r-l 


(A.  7) 


r  I 

Dr(t)  -  -^(tjF^Ct)  +  Dr(t)A*  -  fl2(t)  l  p  LDr(t)  +  H^F’)  " 


n.  (t)  .  Fi(n,(s.  (t))-n,(t))  A'(t-3.(t)) 

n^TtT)  HiG2(F2)  e 


(A.  8) 


then,  for  t>t 
—  x 


/  * 


FI(\ <■)-"!  (tl'v,  Pi<^(*)-''2(t».  A'  (t-a) 


(A. 9) 


L(t)  -  G^D^ (t)  (A.  10) 

where  D. (t  )  is  obtained  from  (A. 7),  (A. 8) ,  and  for  t>t 

X  X  X 

Dj  (t)  -  +  Dj  (t)A'  -  H2<t)Dj+1(t)  +  H[G2(F2)j"lH2  (A*U) 

Dr(t)  -  -fl^tjFjD^t)  +  Dr(t)A'  -  fl2(t)  l  (t)  +  H]G2(F2)r_1H£ 


(A. 12) 
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CM«_bi  As  in  Csss  a,  this  csss  is  tha  saaa  as  Cass  2  up  until  ty 

using  aquations  (2.28)- (2.35)  to  provids  us  with  initial  conditions, 

at  tins  t.  on  B,  and  C. .  For  t,  <t<t 
J  3  3  3  - y 

MM  -  /2<t> ,  W»-V.» 

J  12  G2  ** 

0 

s  (t) 
sJ2(t) 


(A. 13) 


r.(t)  -  Hj^Bj^t)  +  (t)  J  ♦  G^tt)  (A.  14) 

Whsrs  tha  Bj  satisfy  (2.28),  (2.29)  with  sQ(t)  raplacad  by  s2(t).  Nota 

that  B^t^J-O.  Also  tha  Dj  satisfy  (A.7),(a.8)  with  D^t^-o.  iha  only 

naw  aquation  is  for  tha  C.  ovar  tha  tins  intarval  t  <t<t  . 

J  3 - y 

6j  (t)  -  n^tjFjC^t)  +  Cj  (t)A*  -  n2(t)cj+l(t) 

+  7^7t) )  GiGi <r2> 3“X*r2 ^  (#1  (t> (t> )h’sa*  (t’*i (t) > 


n2  (t) 
n2 (*2  (t) ] 


**1  ^1  ^*2  4_i  A' (t-s 

)*  GlGi(Fi)j  H2* 


-•„(t)) 


(A. 15) 
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Cr(t)  -  +  Cr(t)A*  -  f|2  (t)  l  pj_1Cj  (t) 

*l(t)  ,  F2  (l^2  ^S1  (^') }  }  A'  (t-s  (t) ) 

*  ^(tn  gAgi>  •  2  «;•  1 

'A'tl-'iKjItll  ,  A-tt-.2(t» 

n_(s,(t))  a  GlG2(F2)  H2* 


(A.  16) 


For  21 tx*  (t)  and  the  equations  for  and 

are  as  in  Case  3c  -  i.e.  (A.5)-(A.8),  and  C  (t  )-0.  Also,  for  t>t  , 

J  3C  X 

L(t)  -  G^D^t),  and  the  are  calculated  from  (A.ll)  and  (A.12)  . 


Caa*  C:  UP  until  t3  we  are  in  Case  2,  and  thus  use  (2.28)- (2.35) .  over 

the  interval  t3  <t«t  we  have  the  same  situation  as  in  Case  b,  and 

equations  (2.28),  (2.29)  (with  sQ (t)  replaced  by  s2 (t) ) ,  (A.7),  (A.8), 

and  (A. 14)- (A. 16)  apply.  At  tine  t  ,  since  s.  (t  )  -  s„ (t  ),  c, (t  )-0 

C  1C  2  c  j  c 

(see  (A. 13) ) ,  while  (t^)  and  (t^)  take  on  values  calculated  from  the 

previously  mentioned  equations. 

Now  consider  the  time  interval  t  <t<t  .  Note  that  for  t<t 

c - x  c 

Sl(t>>s2(t)'  while  for  fc>tc'  s1(t)<s2 (t).  Thus 
8  (t) 
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0 

r2<t)  .  pi(Vs)"ni(t))  p.(n2(t)-n2(s)) 
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*  yiTftT)  *  »iVa  V 
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E.  (t  )-0 

J  c 


j"l* • ■ ■ it 


(A. 23) 


Dj  (t)  -  -fl^OFjD^t)  +  Dj  (t)A'  -  H2(t)Dj+1(t)  +  H»GJ(FJ)3“1H' 


n,ct)  C  -  \I|-  \«.  \w/  /  —  i  J-  \  W/  an 

•T^Ttir  •  -i-W3  *;• 


f;  (n,  (a.(t))-n,  (t)) 


A*  (t-32 (t) ) 


(A.  24) 


Dr(t)  -  -HjttjF^Ct)  4-  Dr(t)A'  -  fl2  l  +  ^’(F’)*"1^ 


n2(t> 


^ (s2 (t) )  * 


Pi  ">1  (32 <« >  "Hi  <*>  >  H'Gi  <Fi)  r-V.“'  (t“2  M  ’ 

1  1222  (a.25) 


Note  thet  Bj  (t^-O. 

Now  consider  the  interval  t  <t<t  .  In  this  case 

x  y 

s,<t) 

. ^  eel  /m  \  /e\  \  19  I 


r  F.'  (n.  <s)-n,  (t) )  F  (n,  (t)-n,  (s) )  .  . 

J  Gje  1  1  1  H^e  22  2  G2eA  s)ds 


f  G;eFi(ni(S)"ni(t))H*G'eP2(n2<S)'n2(t)>H*eA' (t‘S)ds 


s2(t) 


(A. 26) 


L(t)  -  G£  ID1  (t)  +  E^  (t)  1 


(A. 27) 


where  (t>  satisfies  (A.24) ,  (A.25) ,  with  (t^)  given  by  the  previous 

step.  Similarly,  E  (t  )  is  given  by  the  preceding  step,  but  then  over  the 

3  ™ 


interval  <t<t^ 
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E^tt)  -  -^(tJPjE^t)  +  Ej  (t)A'  ♦  fl2(t)Ej+1(t) 

rL(t)  p:  (n.(«2(t))-n.  (t))  a*  (t-s  <tn 

*  y^TtT)  •  TO  5 2* 


(A.28) 


5 

*r(t>  "  ♦  Er(t)A'  ♦  n2(t)  pj-1Ej 


(t) 


n2  (t) 


^(•2<t))  *  H1H2F2  G2* 


(A.  29) 


Not*  that  E^(ty)-0.  Than,  for  t>ty,  L(t)  -  GJD^t),  whara  (ty)  is 
obtained  from  tha  pracading  aquations,  and  for  t>ty  tha  satisfy 
(A. 11), (A. 12). 


